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ABSTRACT 



nQ I We present calculations of the time evolution of the prompt spectra of gamma- 

[nr^ I ray burst models involving generic internal dissipation regions, including internal 

r^ I shocks, either by itself or in the presence of an external photon source such as a 

7^1 photosphere. The method uses a newly developed time-dependent code involv- 

P I ing synchrotron emission and absorption, inverse Compton scattering and pair 

c/5 I formation. The models reproduce the typically observed Band spectra and their 

generic time evolution, including the appearance of an extra keV-GeV compo- 
nent, whose delay in simple SSC models, however, is only partially able to explain 
K^ I the several seconds observed GeV delays. On the other hand, models involving 

^ I both a photosphere and an internal dissipation region at a larger radius produce 

both an extra GeV component and time delays which are in the range of the 
observations. 



Subject headings: gamma rays burst: general — Radiation mechanisms: non- 
thermal 



1. Introduction 

Most of the prompt emissio n spectra of gamrn a-ray bursts (GRBs) can be described by 



the well-known Band function (IBand et al. Ill993l). which h as a spectral peak around the 



MeV range. Since the EGRET era (JGonzalez et al. 1 120031 ) . distinct spectral components 
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in higher energy ranges have been repo rted. Moreover, brig ht optical emissions during the 
prompt phase, such as in GRB 080319B (JRacusin et al. II2008I ). have also been detected. Such 
lower energy emissions also consist of distinct spectral components. In the popular internal 
shock model of GRB the prompt MeV gamma-rays are attri buted to synchrotron radiatio n 
from electrons accelerated in shocks within the GRB outflow (jPiran II2005I : iMeszaros II2006I ). 
while higher energy photons are attributed to inverse Compton (IC) up scattering of these 



synchrotron photons, a scheme called synchrotron self- Compton or SSC (JMeszaros fc Rees 



19941 : IPilla k Loeb Ill998l : IPanaitescu k Meszaros II2OOOI : iGuetta k Granot 



2003|). Variants 



on this scheme involve having the MeV photons, whether of synchrotron or other origin, 
arising at a different location (usually at a smaller radius) than the locatio n where they 



are upscattered, a scheme refereed to as ex ternal inverse Compton or EIC (JBeloborodov 



! 



20051 : iToma et al. II2009I : iMurase et al.ll201ll ). The generic radiation physics is independent, 
to a large degree, of the specific macroscopic model for the acceleration of the radiating or 
scattering non-thermal particles. Thus, the results we discuss here are not tied specifically 
to internal shock models, but apply also to more generic dissipative models, in which syn- 
chrotron, IC and, in the case of high radiation density also pair formation, are the main 
radiation mechanisms. 

The recent observations over a wide range of photon energies by the Fermi-Large Area 
Telescope (LAT) and GLAST Burst Monitor (GBM) instruments and robotic optical tele- 
scopes such as TORTORA have stimulated research on the temporal evolution of the GRB 
spectra, which is a key for distinguishing various models. Several authors have developed 
time-dependent codes and discuss the spectral evolution in GRBs ( Pe'er k Waxman II2005I: 



Pe'er II2008I : iBelmont et al. II2008I : IVurm k Poutanen II2009I : iBosniak et al. 



20091: 



Daigne et al. 



2OIII ). Some of those studies focus on the effects of the dynamical evolution of the shocked 
region, multiple Compton scattering, or thermal particles etc. An expanding range of phys- 
ical situations needs to be considered, in view of the ongoing observational and theoretical 
progress being made. Here, we improve our previous Mo nte Carlo code, develop ed in a series 
of studies for hadronic and leptonic processes (see, e.g., lAsano k Inoue 1 120071 ). into a time- 
dependent code. Our code is based on the one-zone approximation with the Lagrangian 
specification in the energy space for particles. Although we plan to develop this code to 
deal with the hadronic processes in the future, at this time only the leptonic processes are 
included, as a first step. 

Our goal in this paper is to study the temporal behavior of the radiation field in the 
typical types of geometries generally considered in GRB models. In particular, we study the 
development of the photon spectra as a function of time, as well as the expected light curves 
in specific energy ranges, in order to gain insight into current questions such as the origin 



of the delay between MeV and GeV pulses in bursts observed by the Fermi f lAbdo et al 
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below (JAbdo et al. 



2009a 



2009d ] bl: lAckermann et al.l 20111). or the presence of a dditional spectral components above or 



b|; lAckermann et al.ll201ll ) the usual Band (MeV) spectra known 



from previous Compton and Swift studies. 

In ^ we present the basic physical model used in our numerical, discussing the var- 
ious physical processes involved and the way the time-dependent spectra and light curves 
are computed. The results are presented in ^ In subsection §3.11 we discuss first three 
cases characterized by parameters such as inferred from pre- Fermi GRB observations, in the 
context of an SSC model. These "moderate" models serve as a reference point, providing 
insight into what was missed in the pre- Fermi observations. In subsection §3.2l we explore the 
temporal and spectral properties and predictions for bursts with parameters which might 
characterize some of the more interesting bursts observed by the Fermi LAT instrument, 
again based on an SSC model. In §3.31 we discuss the properties of a model including exter- 
nal photon sources providing a Band function spectrum which is reprocessed and upscattered 
at a different location. We summarize our numerical results in §4] and discuss their relevance 
and implications in ^ 



2. Model and Methods 



In order to simulate the temporal evolution of photon emissions, we d evelop a new 



numerical code by improving our code matured via a series of GRB studies (JAsano 112005 



Asano &: Nagataki I l2006l : lAsano fc Inoue 1 120071 : lAsano et al. I l2009bl ) , in which the photon 
spectrum had been obtained with the steady-state approximation for the photon field. 

Here we consider a shell expanding with the Lorentz factor F from an initial radius 
R = Ro- The calculation of the photon production is carried out in the shell frame (here- 
after, the quantities in this frame are denoted with primed characters). This shell represents 
the internal dissipation region, which in what follows will be referred to as "internal shock" , 
although it could be a generic dissipation region resulting in particle acceleration, e.g., a 
magnetic reconnection zone, etc. Let us adopt notations for particle kinetic energies of 
electrons and photons as ^e = (7c — l)meC^ and e^, respectively. The one-zone approx- 
imation with a constant shell width W = Ro/T is adopted so that the photon density 
n' {e' ) and magnetic field B', etc., in the shell are homogeneous. We refresh the photon 
field with a time step AT' = i?/(100Fc), during which electrons/positrons emit photons 
in the constant photon field. With the same time step we inject electrons intermittently. 
In this paper, we assume a constant injection rate with total energy E^ = TE'^ during a 
timescale t(„5 = W'/c. The injection spectrum is assumed as a cut-off power-law shape. 



mj 



N' ■ -ie' 



cx: £g ^exp (— e^/e^j^^) for e'^ > £^1115 "where e'^^^^ is determined by equating the 



-4- 



cooling timescale t'^^^^ and the acceleration timescale t'^^^ = ^e'^/ceB' (hereafter, the efficient 
acceleration of .^ = 1 is supposed). Using the cooling rates due to synchrotron {£b) and IC 
(£ic), and the heating rate due to synchrotron self absorption (SSA) (^^sa), the cooling time 
is defined as tcooi = ^cK^b + ^ic - 4a)- 

The cooling timescale would be much shorter than the time step AT'. An even shorter 
time step At' = min(AT'/10, t^^Qi/lOO) is used to evolve the electron energy distribution 
during AT' by following each electron with the Lagrangian specification in the energy space. 
In every time step At', we accumulate emitted photons, and add them to the photon field 
n' [e' ) in every AT'. In the following subsection we explain the detailed methods to estimate 
t[,QQi and simulate the emission processes. The methods are essentially the same as our old 
code. 



2.1. Physical processes 



2.1.1. Synchrotron emission 
The energy loss rate due to synchrotron/cyclotron emission is 

£b = -—dT^cB'^-ilPl sin^ a, 

47r 



where o"t is the Thomson cross section, and a is the pitch angle (JRybicki fc Lightman Ill979l ). 
In the ultra-relativistic limit (7e 3> 1), the energy spectrum of synchrotron emission per unit 
time per photon energy is written as 



PB{£e,£-y) 



VSe^B 



sma 



X 



2TThmeC^ 

2meCe^ 
3'y'^ heB sin a' 



Fix) 



(2) 
(3) 



where F{x) is the synchrotr on function. In the noii -relativistic energy range, instead of 
the exact formula (see , e.g., iPe'er fc Waxman 1 12005| ). we adopt a simple approximation 



flGhisellini et al. lllQQSf ) for /3e = ^/l - \hl < 0.89 as 

mecVT^ 72/32 [2(1 - e^/es) 



PB{Bc,e.y) 

£b 



Sirhe 1 
heB 



37e'A 
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(4) 
(5) 



for the isotropic distribution of electrons. 
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2.1.2. Inverse Compton 



In order to calculate emission due to photon scattering with electrons, we prepared 
tables in advance. Since we plan to develop this code to simulate optically thick plasmas 
in the future, the tables cover a wide energy range from the Thomson regime to the Klein- 
Nishina limit. In the laboratory frame, the photon distribution is assumed to be isotropic. 
Integrating over incident and scattering angles with the Monte Carlo method, we obtained 
the tables of the emission spectrum as below. Here, we denote the quantities in the electron 
rest frame by characters with tilde. Given the energies of an electron e^ and an incident 
photon e^ in and a cosine of incident angle /lin, the Lorentz transformation yields the photon 
energy in the electron rest frame e^^^. The energy of the scattered photon depends on the 
cosine of the scattering angle fi as 



'7, out 



'7, in 



1 + {e-^.\n/meC^){l -/i)' 
The scattering probability against jl is 



V-^(x 



'7, out 



-7,in 



2 r ~ 

e 



7, in 
^7, out 



-7, out 



-7, in 



;i-/i^) 



(6) 



(7) 



which should be normalized by the total cross section as 



C"KN 



-CTT 

ln(l 



1 + a; f2x{l + x) 



2x) 



l + 2x 
l + Sx 



ln(l + 2x)) + 



X 



(9) 



2x 

~ I 2 

(JRybicki fc Lightman Ill979l ). Reverting to the laboratory frame again, we obtain the tables 
of the emission spectrum normalized by the photon density in a form of 

-< ic(^c) ^7, in; ^7,out)/''^7(£^7,in)- (10) 

At the same time we obtain the tables of the "photon-absorption" probability for incident 
photons, which assures the conservation of the photon number in the scattering process. The 
Monte Carlo integral over the angles brings us another table of the average value. 



Ae^dKN 



Pe/^inJ 1,^7, out 



-7, in 



)0"KN, 



:iii 



which is useful to calculate the energy loss rate as 



8 



ic 



c / de^n^{e~^)^e^a}i^{ef,,e^) 



(XC'vCZc 'yout-'ICv^e; ^7; ^ 



7? t^^outj 



(12) 
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2.1.3. Synchrotron self- ah sorption 



Adopting; the synchrotron/cyc lotron formulae in §2.1.H the cross section of SSA (JRybicki fc Lightman 



19791 : iGhiselhni fc Svensson Ill99l[ ) is written as 

(7abs - l)mcC^ 



C^SA 
tabs 






PsK^ci^'y, 



(13) 
(14) 



A fraction 1 — exp (— cAt' j decUcasA) of photons will be absorbed every time step. The 
accumulated amount of absorbed photons is counted in every AT'. The heating rate of 
electrons due to SSA is 



^SA = / dSjCSyasAn^i^ 






(15) 



2.1.4. Electron-Positron pair creation 

In order to take into account 77-pair creation, every time step AT' we eliminate a 
fraction of photons, 1 — exp {—t±), with 



T"±l^7,lJ 



cAT' / dej^inil - /ii: 



"^71^7,2 

47r 



-o-±, 



where 



<y± 



16 



(Jt(1 -y 



(3-y^)lnl±^-2|/(2-y2) 



y = 1- (2m,c )/[e^,ie^,2(l-/Um)] 



(16) 

(17) 
fl8) 



( iBerestetskii et al.l Il982l ) . The integral over the photon energy and incident solid angle 
(ifiin = 27r(i/iin is numerically calculated with the isotropic photon distribution. We simul- 
taneously inject secondary electron/positron pairs. To save computational cost, we adopt a 
simple approximation for the energy of the secondary pairs as e^ = (£7,1 -1-6^^2)/'^ — fn^c^ and 
neglect the effect of electron-positron pair annihilation. 



2.1.5. Shell expansion 

We take into account the effect of adiabatic cooling, though it may be a minor effect in 
GRB physics. The adiabatic cooling of the collisionless plasma is not trivial. The motion 
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of electrons/positrons is controlled by only magnetic fields. If a homogeneous magnetic field 
evolves according to volume expansion, only the momentum component perpendicular to the 
magnetic field would be affected. However, the magnetic field in the GRB sources may be 
highly entangled so that the pitch angle diffusion cannot be neglected. Thus, in this paper, 
we use a simple isotropic approximation for the "momentum", \Je'^ + 2mf.(?E'^c~^ oc ^'~^/^, 
which is applicable from non-relativistic (£g oc l^'~^/^) to ultra-relativistic (£g oc y~^/^) 
limits. The constant shell width W assumed in this paper implies that the volume evolves 
as V = AirR^W oc R^. The volume expansion may imply evolution of the magnetic field, 
but we assume constant B' for simplicity in this paper. 



2.1.6. Photon escape 

The photon density n'^ refreshed every time step gives us the photon escape rate per 
unit surface as n'^c/A. The photons escape from both the foreside and backside surfaces of 
the shell (the total surface: 2 x AttR^). Every time step we extract a fraction of photons, 
2nR^cAT'/V' = cAT'/{2W'). 



2.2. Light Curve 

We consider photons escaping from the shell with latitudes 6 < 6jct to take into account 
the jetted structure, but the maximum latitude is taken as large as 6'jet = 5/r so that the 
effect of collimation does not affect observation very much. The energy and escape time t' of 
photons coming from a latitude 6 are transformed into energy- and time-bins in observer's 
frame with 

e' 

and 

Us = {l + z) [(1 - /3shCose)n' + Roil - cose)/c] , (20) 

respectively, where /3sh = -\/l — l/F^. Every time step AT' the spectral number of escaped 
photons per unit solid angle and surface is estimated with the photon density n'^ as 



where we count photons escaped from both foreside (cos 6*' > 0) and backside (cos 6'' < 0) of 
the shell with the Lorentz transformation for 9 

cose-/3sh i^^. 

1 — Psh COS 6 

For high latitude emissions with cos 9' < 0, the shell thickness adds an extra "time delay" 

WT 
Aibk = (1 + ^) ^^ (cos 9 - /3,h) , (23) 

in Eq. f l20|) . The number of photons escaping from a surface element A5" 

dN'ie') , 

-^A5', (24) 

is transformed to the frame of the central engine with 

dS' = dS = 2nR^ sin 9d9, dVt' = -r- ^ — , dN' = dN. (25) 

r^(l-/3shCos6')^ 

In each corresponding time-bin, the observer counts the contribution from AS* for the photon 
number per unit surface as 

~ds~^, ^ Dldms ' ^^^^ 

where the proper distance Dp = -Dl/(1 + z) is used, because we have already included the 
effects of the energy shift and time-dilation due to cosmological redshift in eqs. (TT9|) and 
( 120|) . Integrating over the surface with < ^ < ^jct, we obtain the total fluence for each 
time-bin. 



As was calculated in iMurase et al.l ( 120071 ). the photon absorption due to extra galactic 



background light (EBL) is estimated by accumulating the optical depth between the source 



and observer with the EBL model of iKneiske et al. I (12004 ) (best-fit) 



3. Results 

The model parameters are five: the initial radius Rq, the bulk Lorentz factor F, the 
comoving magnetic field B', the electron injection energy E^,, and the minimum electron 
energy e'^jn- While we can consider more complicated situations by evolving e.g., B', W, 
F, the electron injection rate, etc., with R, we adopt here for simplicity the constant values 
of the parameters as mentioned in [|2j The other fixed parameters are the jet opening angle 
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^jct = 5/r, electron injection index p = 2.5, and the parameter for the electron acceleration 
timescale ^ = 1. 

Since we make no assumption about the energy density of the "thermal" components, 
the conventional parameters of th e energy fraction to the total internal energy density, such 



as ee and e^ (see the definitions in lMeszaros II2006I ). are not specified. But assuming ee = 0.5 



with the electron energy density Ec/{4:7tRI), hereafter we will take e^ as a referential value 
for each parameter set. The Thomson optical depth is another criterial parameter. The 
electron number density in our simulations can be approximated as 

n' = (^-^)^- (27) 

In all our simulations, the Thomson optical depth tt = n'eO"T-Ro/r is much smaller than 
unity. 



3.1. GRBs with "Moderate" Parameters 

In this subsection we discuss three cases labeled Runl, Run2, and Run3, whose param- 
eters are indicated in the following figures [TJIHI In Figure [H we show the evolution of the 
spectral photon density in the shell frame for Runl. The model parameter set is one to be 
considered as fiducial values before the Fermi satellite was launched (if e^ = 0.5, e^ = 0.0025; 
Tt = 0.03). As the electron injection proceeds from R/Rq = 1, the photon energy density 
gradually grows. In this case the energy density of the magnetic field is ~ 4 x 10^ erg cm~^, 
which can be compared with the vertical axis in Figured] At an early stage (see R/Ro = 1.1 
in Figure [1]) synchrotron emission is the dominant cooling process for injected electrons so 
that the spectral peak associated with e^j^ is seen at a few keV. The relatively low photon 
density at this time leads to a higher cut-off energy (~ 10 MeV) due to 77-absorption. As 
the photon energy density increases, the IC component gradually becomes prominent (see 
R/Rq = 1.2 and 1.5), and the cut-off energy due to 77-absorption decre ases. This late 



growt h of the IC component has been observed also in the simulations of iBosnjak et al. 



(120091 ). The IC component makes the spectrum so hard that the synchrotron peak gradually 



disappears, and instead, 77-absorption produces a spectral peak at a few hundreds keV. 

Below the synchrotron peak energy at a few keV, the synchrotron photons from cooled 
electrons of e'^ < e^j^ make a power-law spectrum with the photon index ~ —1.5. The 
l ow-energy cut-off due to SSA is clearly seen at ~ 0.1 eV. The electron heating via SSA 



(JGhisellini et al. Ill988[ ) makes a spectral bump just above the cut-off energy. At the end of 



the electron injection {R/Rq = 2), the superposition of the synchrotron and IC components 
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shows a hard power-law spectrum above the synchrotron peak energy. Considering the 
Klein-Nishina effect, the majority of seed photons for the IC process would be low-energy 
photons at ~ 1 eV. Above the cut-off energy at a few hundreds keV, photon production 
and 77-absorption balance with each other, and consequently generate a steep power-law 
spectrum. After the electron injection stops, the shell expansion and photon escape dilute 
the photon density (see the line of R/Rq = 2.5). In this stage, the photon production is 
practically terminated so that our calculation based on the one-zone approximation shows 
a sharp cut-off above the 77 cut-off energy. The photon energy density in the shell frame 
at R/Ro = 2 is 1.4 X 10^ erg cm""^, which is close to the simple estimate (note that W' is 
constant) EJ{A'kR^Rq) with R = 2Ro, 2.0 x 10^ erg cm'^. 
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Fig. 1. — Temporal evolution of spectral energy density of photons in the shell frame for 
Runl. The model parameters are denoted inside the figure. The electron injection is ended 
at R/Ro = 2.0 (thick line). 



According to the method described in §2.2[ we plot light curves based on the spectral 
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Fig. 2. — Temporal evolution of the observable spectral photon flux for the model described 
in Figure [1] (Runl). The time-bins are the same as in Figure [3l and the denoted times are 
the center for each time-bin. 



flux (the unit is erg cm~^ s~^ eV~^) assuming the source redshift of 2; = 1 in Figure |3l The 
time-bins are logarithmically divided (0.25 dex normalized by (1 + z)RoT~'^c~^), and the 
flux is averaged over each time-bin. In the flgure, the GeV, 100 MeV, and 100 keV light 
curves are plotted. The keV light curve practically overlaps with the 100 keV lightcurve, 
and th e light curve shapes seem to reproduce well the typical FRED (fast rise, exponential 
decay; iFenimore et al.lll996l ) shape, which is frequently found in many GRBs. Reflecting the 
temporal decrease of the 77 cut-off energy in the source, the GeV flux peaks earlier than the 
lower energy fluxes. This corresponds to what one would get from an extrapolation o f the 
"positive spectral lags" , observed in many of the pre- Fermi bursts (JNorris et al.l Il996l ) (see 
also §3.2p . to the GeV energy range. This behavior is explicitly interpreted by the evolution 
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of the photon spectrum shown in Figure [2j Before the peak time of the 1 keV-100 MeV 
fluxes (~ 58 ms), the spectrum seems to reflect the evolution of the spectrum in the source 
frame. During the decaying phase, the spectrum gradually shifts to lower energy owing to 
the emissions from higher latitudes of the shell (curvature effect). We can see that GeV 
emission ceases faster than lower energy bands, which explains the narrow pulse shape in 
the GeV light curve. As shown in Figure HI however, the GeV fluence is dim compared to 
the MeV-fluence in this parameter set. 

Normalized Flux 

I I I I • • • ' 



z=l - 
Ro=10''^cm, r=300, B'=10'^G 
E,=10^' erg, eU= 1.26 GeV 




0.4 t[s] 



Fig. 3. — Light curves for the model described in Figure [T] (Runl). The plotted values are 
normalized by the maximum flux. We smoothly join the fluxes averaged over respective 
time-bins with solid curves for reference. 

We have shown the results for Runl as an example that shows a remarkably suggestive 
evolution of the spectral shape. However, the high-energy spectral index is harder than 
the typical observed ones. Hence, we consider two different cases, in which the internal 77- 
absorption for GeV photons is not efficient, and the IC component is not so prominent; Run2 



-13- 



'1 

8F(8) [erg/cm ] 

rnnif 




Fig. 4. — Summarized observable fluences for Runl-Run3 with z = 1. The dashed hnes 
denote the cases neglecting 77-absorption due to EBL. 



with large T and Rq (Figures [5] and (6] if eg = 0.5, e^ = 0.06; tt = 3 x 10~® ), and Run3 with 
moderate F and large magnetic field (Figures [7] and [HI if Ce = 0.5, e^ = 0.025; tt = 9 x 10"^ 
). The relatively stronger magnetic field than that in Runl and the Klein-Nishina effect 
make synchrotron radiation the dominant cooling process even though ec > e^. In these 
cases the evolution of the fluxes (Figures [5] and [7]) are monotonic compared to Runl. 

In Figure [5] the low-energy spectrum below epcak — 300-400 keV becomes soft with time. 
The slow electron cooling due to the low magnetic field in Run2 causes the slow evolution of 
the low-energy spectrum. The typical energy of electrons emitting 1 eV synchrotron photons 
is ~ 8 MeV. The synchrotron cooling time for such electrons ~ 550 s is slightly longer than 
t^Qj ~ 100 s, although the IC cooling is marginally important even after t' = t[^y This is the 
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reason for the slow spectral evolution, but note that t'^^^^ for electrons of e^j^ is much shorter 
than the dynamical timescale. So almost all the energy injected in the shell is released as 
photons even in this case. For Run2 the SSA frequency is well below the optical band so that 
we can expect bright optical pulses as shown in Figure HI While the light curves for keV-100 
GeV regions are almost the same as the 100 keV light curve in Figure [6], the optical (~ 1 
eV) light curve shows a significant delay relative to the higher energy bands. This delay is a 
direct consequence of the gradual softening of the low-energy spectrum. The broader pulse 
profile of the 1 eV light curve compared to that of the 100 keV light curve seen in Figure 
is due to the emission from higher latitudes. As shown in Figure HI the intrinsic spectrum 
for Run2 extends as far as the 1 TeV range, where the IC component dominates, but such 
photons are absorbed by EBL in our assumption of 2; = 1. 
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Fig. 5. — Temporal evolution of the observable spectral photon flux for the model Run2 (see 
also Figure EJ. The time-bins are the same as in Figure El and the denoted times are the 
center for each time-bin. 
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Fig. 6. — The light curves for the model Run2, whose parameters are denoted inside the 
figure. The plotted values are normalized by the maximum flux. We smoothly join the fluxes 
averaged over respective time-bins with solid curves for reference. 



The time-integrated spectrum for the more compact case of Run3 is similar to that of 
Run2 (see Figure H]) so that we can expect an optical pulse and 100 GeV photons possibly 
detectable with Cerenkov detectors such as MAGIC or CTA for those two cases. The high- 
energy cut-off in the spectrum is mainly due to the internal 77-absorption rather than due 
to the EBL, unlike for the case of Run2. The relatively fast electron cooling due to the 
high magnetic fleld in Run3 promptly generates a soft photon spectrum with a photon index 
~ —1.5 below ^pcak- Thus, in this case no signiflcant softening is found during the pulse rise 
phase (< 0.58s) unlikely Run2. On the other hand, the smaller radius than that in the case 
of Run2 yields a slight evolution of the high energy cut-off due to 77-absorption. Therefore, 
the peak of the 100 GeV light curve appears earlier than that of the lower energy bands, 
while the light curves for keV-GeV regimes almost overlap with the 100 keV lightcurve. The 
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Fig. 7. — Temporal evolution of the observable spectral photon flux for the model Run3 (see 
also Figure IH]) • The time-bins are the same as in Figure El and the denoted times are the 
center for each time-bin. 



heating of electrons due to SSA (JGhisellini et al. Ill988l ) is so effective that a slight bump 
appears at a few eV in the spectrum in Figure HI The high-energy bump, prominent above 
100 MeV comp ared to Run2, rnay be due to secondary electron-positron pairs produced via 
77- absorption (JAsano fc Inoue 1120071 ). 



Above the 77 cut-off energy, the intrinsic spectral shape for Run3 is well approximated 
by a power-law rather than by an exponential cut-off. This steep power-law shape is also 
seen in the spectrum of Runl. If we take into account the spatial structure o f the source, a 



smoothly broken po wer-law would be expected as a 77-absorption signature ( JBaring 112006 



Granot et al. 1120081 ). We see that, even with our one-zone approximation, the temporal evo- 
lution of the spectrum generates a broken power-law shape in the time-integrated spectrum. 
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Fig. 8. — Light curves for the model Run3, whose parameters are denoted inside the figure. 
The plotted values are normalized by the maximum flux. We smoothly join the fluxes 
averaged over respective time-bins with solid curves for reference. 



The poor statistics in the spectral data of GRB 090926A ( jAckermann et al.ll201lh . for which 
the spectral break is reported around 1.4 GeV, makes it difficult to distinguish whether the 
spectral shape is one of a broken power-law or an exponential cut-off. Thus, future detections 
with better statistics are needed of bright bursts with Fermi or with Cerenkov telescopes, in 
order to determine more precisely the spectral shape above the cut-off energy. 

Our results show that the curvature effect yields a spectral lag between the optical 
and the keV-MeV bands. In the keV-MeV energy range, several authors have reported time 
delays in th e arrival time of lower energy photons relatiye to higher energy photons ("positive 



lags" , e.g.. iNorris et al. 



1996 



20051 : iHakkila et al. II2008I : lArimoto et al. Il2010[ ). These spectral delays may be explained by 



Wu fc Fenimore I l2000l : iHakkila fc Giblin I l2004l : IChen et al. 
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the cu rvature effect fjloka fc Nalcamura ll200ll : IShen et al.ll2005l : iLu et al.ll2006l ). As lLu et al. 
( 120061 ) discussed, the curvature effect predicts that a smaller T tends to yield a larger lag, 
i.e., softer photons arrive increasingly later. The optical delay in our results is consistent 
with this tendency, the smaller F in Run3 resulting in a broader 1 eV light curve than that 
in Run2 (note that the optical delay in Run2 is mainly due to the slower electron cooling). 
However, the observed energy dependence of the l ag does not always fol low the prediction 
of the simple curvature effect ( IZhang et al.l 120071 : lArimoto et al. 1120101). We should also 



take into account the tempo r al evolution of the spectra l shape (JKocevski fc Liang 1 12003 



Daigne fc Mochkovitch Ill998l . l2003l : iBosniak et al 



20091 ). It would be interesting to fit the 



observed spectral lags by evolving some of the model parameters, such as the magnetic field, 
although it is not our purpose to discuss such effects in this paper. 

The broad pulse profile in the optical band means that the optical pulses tend to overlap 
each other. The observed broader pulses and lo nger tail in the optica l light curve compared 



to the gamma-ray light curve in GRB 080319B (JRacusin et al. II2008I ) is consistent with this 



tendency, although an extra spectral component is required to reproduce the optical fiux in 
that burst. Recalling that the slow electron cooling is a key to the optical lag in Run2, we 
emphasize that the spectral evolution should be investigated not only in the GeV band but 
also in the optical band in order to probe the prompt emission mechanism. 



3.2. "Fermi-LAT" Case 



The Fermi satellite dete cted GeV photons f rom several very bright bursts (-Ejsn > 10^^ 
erg) such as G RB 080916C JAbdo et al. I l2009ch. GRB 090902B jAbdo et al. Il2009bh . and 
GRB 090926A (jAckermann et al.ll201ll ). In such bursts the onset of the GeV emission is 
delayed with respect to the MeV emission, i.e., they have (in this energy range) negative 
lags. In addition, some of them also have an extra spectral component above a few GeV, 
distinct a nd additional to t he usual Band function which dominates in the MeV energy 
range. As ICorsi et al.l ( 120101 ) discussed, the GeV emission may be due to IC emission from 
internal dissipation regions. In this subsection we choose a set of parameters, which may 
be relevant for such bright GRBs accompanied by an additional GeV emission component, 
although we do not attempt to perform fits to the spectrum of any specific GRB. 

In Figure [HI we plot the evolution of the spectral photon density in the shell frame for 
Run4, in which we have adopted a very large injection energy [Ec = 10^^ erg) and a high 
Lorentz factor (F = 1000). If we adopt eg = 0.5, B = 100 G corresponds to e^ = 5.4 x 10""^. 
The Thomson optical depth is low enough, tt = 3 x 10^^, even for the large Eg. 
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Fig. 9. — Temporal evolution of the spectral energy density of photons in the shell frame for 
Run4. The model parameters are denoted inside the figure. The electron injection is ended 
at RjR^ = 2.0 (thick line). 



Initially the synchrotron component at ~ keV dominates the spectrum, but the IC 
component grows as electrons are continuously injected. In the later stage the injected 
electrons cool mainly via IC rather than synchrotron, so the synchrotron component starts 
to decay earlier than the IC component. At the end of the electron injection [RJRq = 2) the 
cooled electrons are still relativistic, because the energy density of the magnetic field is low 
(~ 400 erg cm~^). Equating the synchrotron cooling timescale and the dynamical timescale 
W'/c, the energy of the cooled electrons is 
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The contribution of IC cooling with the Klein-Nishina effect can also be comparable to the 
above estimate, as the obtained spectra show; the photon energy densities of the IC and 
synchrotron components are comparable. Actually, our numerical results indicate a spectral 
bump in the electron distribution around 100 MeV at R/Rq = 2. Even after the end of 
electron injection, such cooled electrons continue emitting synchrotron photons. The cooling 
due to IC gradually becomes inefficient as the seed photon density around keV energies 
decreases. The typical synchrotron photon energy from the cooled electrons is 



'7, cool 



SnheB' 
SrrieC 

0.2eV 
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In the spectra for R/Ro = 2 and 2.5 we can see a spectral bump below 
attributed to the late synchrotron emission from the cooled electrons. 



(29) 



1 eV, which is 



Assuming a redshift z = 4.35 such as that of as GRB 080916C, we show the evolution 
of the spectrum in the observer frame for this Run4 case in Figure dUi The IC component in 
the GeV-10 GeV range grows with time, and the late synchrotron emission discussed above 
produces a spectral break at keV ~ Te' ^^^^ The spectrum seems to evolve from the initial 
Band-function to a power-law dominated shape in th e later stages. This spectral evolution 



is int erestingly similar to that seen in GRB 090510 ( lAbdo et al. Il2009al : lAckermann et al. 



2010l ). although that burst is a short GRB. The spectrum evolves to a power-law-like shape 
well after the peak time of the pulse (typically t = 2.6 s in Figure [101 see also Figure 
[T3!) . If many pulses overlap each other, a dominant pulse may prevent the detection of the 
underlying power-law-like component. Therefore, it is interesting that GRB 090510 is one 
of the short GRBs, which have only one or two pulses in most cases. In many Fermi-LAT 
GRBs, the spectra tend to evolve toward a simple power-law spectrum extending as far as 
GeV in the later stage of the prompt emissions. Even for such long GRBs, the spectral 
evolution seen in our simulation would be applicable, because active pulses may gradually 
disappear in t he later stage. While such spectra may be explained by the early onset of 
the afterglow fJGhisellini et al. Il2010l : iKumar fc Barniol Duran Il2010[ ). a combination of the 
late synchrotron emission and a decaying IC component can also produce a power-law-like 
spectrum. 

The time-integrated spectrum in Figure [11] shows extra components below keV and 
above 100 MeV , comparable to those seen in GRB 090510 (lAckermann et al.l l2010l ) and 
GRB 090926A (lAckermann et al.ll201ll ). To explain the observed extra components that 
dominates the emission both below 10 keV and above 100 MeV, a single contribution from 
an IC component is not sufficient. However, the late synchrotron emission can assist in 
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Fig. 10. — Temporal evolution of the observable spectral photon flux for the model Run4 
with z = 4.35. The time-bins are the same as the logarithmic time-bins in Figure HS} and 
the denoted times are the center for each time-bin. 



extending the extra component into the lower energy range, as shown in Figure [TTl It 
seems that a single power-law component with a photon index —1.82, which is close to the 
index —1.79 in GRB 090926A, overlaps the usual Band function that dominates the emission 
around MeV, despite the fact that the two spectral excesses in the low and the high energy 
regions have different physical origins. 



In Figure [12] we plot the light curves based on the spectral flux with linear time-bins, 
in order to distinguish the peak times. The light curves show a delayed onset of the optical 
emission, which is partially due to the late synchrotron emission. The late growth of the 
IC emission makes th e growth of the early GeV light curve slower compared to the 100 
keV emission (see also lBosnjak et al. II2009I ). the GeV peak time appearing delayed ~ 0.4 s 
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Fig. 11. — Time-integrated observable flux for Run4 with z = 4.35. The dashed hnes denote 
the cases neglecting 77- absorption due to EBL. 



relative to the 100 keV peak time. One needs however to take into account the low photon 
statistics above 100 MeV in the GRB observations with Fermi. Thus, in Figure [13] we plot 
the expected photon-count rates integrating over the 100 MeV - 1 GeV and the > 1 GeV 
ranges. Here we have simply assumed constant effective areas of 3,000 cm^ and 7,000 cm^ for 
energies below and above 1 GeV, which mimic the effective are of the Fermi-LAT detector 
( jRando et al. 1 12009| ) . While the count rates with logarithmic and linear binning give us 
rather different impressions, the flgure shows that we expect from this simple one-zone SSC 
model a detection of one photon with the Fermi-LAT, delayed at least ~ 1 s after the onset 
of the 100 keV or MeV emission. On the other hand, the observed delay of the GeV emission 
seen in GRB 080916C is ~ 4 s, which is substantially longer than 1 s. Although the GeV 
delay due to the growth of the IC component would be within the approximate timescale of 
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Fig. 12. — Light curves for the model Run4 with z = 4.35 based on the spectral flux. The 
time-bins are linearly divided. The plotted values are normalized by the maximum flux. We 
smoothly join the fluxes averaged over respective time-bins with solid curves for reference. 

the keV-MeV pulse width, this delay appears able to explain only a fraction of the observed 
delayed onsets of the GeV emission in long bursts. 



3.3. External Photons 



As discussed in the last subsection, the IC emission tends to be delayed relative to the 
main keV-MeV emission, but the lo nger delay than the pulse timescale of the MeV emission 
such as observed in G RB 080916C (jAbdo et al. I 2009d ) is not explained by this effect only. 
As a possible solution, iToma et al. I (J2009|, 120101 ) pointed out that the delayed GeV onset 
can be reproduced by the introduction of external photons such as the X-ray emission from 
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Fig. 13. — Expected photon counts per each time-bin above 100 MeV and GeV for the 
model Run4 with z = 4.35. The thick (thin) histograms correspond to logarithmic (linear) 
time-bins. For comparison we also plot the lOOkeV- light curve (arbitrary unit), whose shape 
is the same as that in Fig. [T^i 



the cocoon (see also iLi Il2010l : iMurase et al.ll201ll ). The geometrical configuration in those 
models, namely the spatial separation between the source of the external photons and the 
site of the internal shock or generic dissipation region can cause the dela yed arrival of the 



up-scattered external p hotons. In particular, the photospheric emission (jPaczynski I Il986 



Meszaros fc Rees Il2000l ) can provide the external seed photons for the IC emission from elec- 
trons accelerated in regions outside t he photosphere. As suggested, e.g., for GRB 090902B 
(JRyde et al. Il2010l : IZhang et al.ll201ll ). many authors have discussed the photospheric emis- 
sion as the main MeV component. 



In this subsection we assume that a quasi-steady emission provides the main MeV 
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component, emanating from an inner region, whi le an outside dissipa tion region producing 
accelerated particles upscatters these photons, as iToma et al. I ( 120 lOl ) assumed. The quasi- 
steady emission may originate from the photospheric emission, but we do not specify its origin 
here. The mod el details are as fol lows: The spectral shape of the quasi-steady emission is a 
Band function (JBand et al. Ill993l ). whose parameters are the luminosity Lseed, peak energy 
epeak, and spectral indices a and /3. Here, the external emission is assumed to start at a 
time ~ -Rq/c before the electron injection. In this Band photon field, accelerated electrons 
are injected in a shell expanding from R = R^iw the same way as we did in the previous 
subsections. In the shell frame, the external photon field should be anisotropic, because the 
Band component may be emitted (depending on the specific model) from a radius ^ -Rq; 
and its source may have a different bulk Lorentz factor than the F of the shell. Although our 
one- zone numerical code cannot fully include the effects of such an anisotropy, we endeavor 
to take it into account partially as follows. As long as the electron distribution in the shell 
frame is isotropic, the total IC emissivity, even for anisotropic seed photons, is written with 
the average intensity of the seed photons as 



J{e^) = 



Air 



(30) 



(JAharonian fc Atoyan Ill98ll : iBrunetti 1 120001 ) . where Useed{£-y) is the spectral energy density 
of the seed photons. However, reflecting the distribution of the seed photons, the resultant 
emissivity also becomes anisotropic. If the seed photons are emitted from a radius <^ Rq, 
we can consider highly beamed seed photons that radially propagate even in the shell frame. 
In this case the anisotropy in the IC emissivity may be approxi mately included by adding 



an extra factor (1 — cos 6') in Eg . (12T|) in the Thomson limit (jWang &: Meszaros 1 12006 



Fan et al.l l2008l : iToma et al. 1 12009| ) . This is because electrons moving in the same direction 
as the photon beam have a smaller probability to interact with those photons, while con- 
versely the head-on collisions between electrons and photons is much more efficient. In this 
approximation the IC emission becomes maximum from the marginally high-latitude surface 
of the shell, rather than from the on-axis surface. For a totally beamed photon field, its 
spectral shape in the shell frame is obtained by shifting the above Band function to a lower 
energy by a factor of (1 + /3sh)F. The average intensity of the Band component in the shell 
frame is calculated with a normalization of the energy density 

for the totally beamed case. This provides an extra factor (1 + /3sh), compared to the isotropic 
seed-photon case. 

The results shown below are obtained with the fully-beamed approximation mentioned 
above. While expedient and useful, some of the weaknesses of such an approach are listed for 
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Fig. 14. — Temporal evolution of the observable spectral photon flux for the model Run5 
including the external Band component with z = 4.35. The model parameters are denoted 
inside the flgure. The fluxes are averaged over the logarithmic time-bins in Figure [15] and 
the denoted times are the center for each time-bin. The solid lines are calculated based 
on the fully-beamed approximation, while the thin dashed line represents the results in the 
isotropic approximation. 



completeness below. The method in our code corrects for anisotropic effects all of the emis- 
sions from the shell, including synchrotron, which should be independent of the anisotropy 
of the extra photon fleld. Moreover, the seed photons for the IC include not only the exter- 
nal Band component but also the photons emitted in the shell itself (SSC or secondary IC 
emissions). Such secondary seed photons may soften the IC anisotropy. The Klein-Nishina 
effects may also result in deviations from a simple anisotropic factor (1 — cos^'). Even if the 
anisotropy of the emissivity in the shell is simply expressed by this factor, the 6''-dependence 
of the escape timescale should yield the temporal evolution of the anisotropy of photon distri- 
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bution. Photons with small | cos6''| remain in the shell longer than photons with | cos6''| ~ 1. 
Thus, the angular distribution should evolve even for the constant angular dependence of the 
emissivity. Nevertheless, the one-zone and fully-beamed approximations force the angular 
distribution to remain as (1 — cos 6'). Given the average photon energy e:^ in the shell, our 
treatment results in the average momentum of the escaped photons to become —e^/{2c), 
although that of the photons in the shell is — e^/(3c). Owing to the negative average mo- 
mentum, the fluence for the observer decreases by a factor two compared to the isotropic 
approximation. 



The anisotropy of the Band component may affect 77-absorption as well f lZou et al. 



I2 



201 ll : IZhao et al.ll201ll ). The high-energy cut-off can be higher than that calculated here. 
In this case too, however, the target photons against 77 are not only the external photons. 
Since in any case our one- zone approximation is not suited for including the effect of the 
anisotropy in the absorption, we do not include these effects here. 

We adopt the parameters of the Band spectrum as Lsccd = 10^^ erg s~^ and epcak — 1 
MeV in the observer frame, a = —0.6, and /3 = —2.6 in the simulations of Run5 (Figures [T^ 
and fTH]) . The parameters for the internal shock are similar to those in Run4 except for a small 
e^jn (if ee = 0.5, €3 = 0.0018; tt = 2 x 10^^). The small e'^^^ leads to a large implied proton 
energy. The proton density n'p would be the same as the electron density estimated from Eq. 
f l27|) . Even if we neglect the thermal energy, the proton energy Ep ~ TmpC^n'pV ~ 1.9 x 10^^ 
erg is fairly large. In Figure [H] we show the spectral evolution for the observer. We see that 
the synchrotron and IC emission in the shell contributes to a wide energy range, while the 
Band component is prominent mainly in the MeV range. 

In Figure [15] we show the corresponding light curves. We see that the GeV and eV 
emissions originating from the (outer) shell are delayed ~ (1 + z)Rq/T'^ /c ~ 1 s relative to 
the MeV emission, which is "external", i.e., it originates from a region which is at smaller 
radii than the shell. We should note that the start time of this external emission can be 
arbitrary adjusted, depending on the model, so that the delay timescale could be longer than 
shown by our results. Since we assume a quasi-steady Band component, the MeV emission 
continues even after the end of the GeV emission. For comparison, we have also plotted 
in Figure (TS] the light curves using the isotropic approximation for the photon emissions 
(thin lines). Since the marginally high-latitude emission contributes the most to the flux 
in the fully-beamed approximation, the light curves show larger delayed onsets and longer 
tails relative to those in the isotropic cases (~ 0.3 s). The Lorentz boost of the anisotropic 
emissions reduces the fluence to half of that of the isotropic emission owing to the negative 
average photon momentum in the shell frame. Thus, those lower fiuences and longer emission 
timescales in the fully-beamed approximation result in the lower fluxes shown in Figure [1] 
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Fig. 15. — Light curves for the model Run5 including the external Band component with 
z = 4.35 based on the spectral flux. The normalizations are arbitrarily shifted to un-merge 
the lines. The solid lines are plotted based on time-bins linearly divided (0.0423 s). The 
fluxes averaged over logarithmic time-bins are also plotted. The fully-beamed and isotropic 
approximations are given by the thick and thin lines, respectively. The peak fluxes for the 
thick and thin lines are aligned in the GeV and eV light curves, while the MeV light curves 
are normalized with the external photon flux. 



The origin of the eV emission is basically synchrotron, so the light curve in the isotropic 
approximation (thin line) may be more appropriate for this. Even for the IC emission, the 
actual flux evolution may be between the cases represented by the fully-beamed and the 
simple isotropic approximation. 

The internal shock produces the synchrotron and IC components at ~ 10 eV and ~ 10 
GeV, respectively. The sync hrotron peak appears able to explain the optical extra component 



observed in ORB 080319B (iRacusin et al. II2008I ). In the later stages, the seed photons for 
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the IC process are not only the Band component but also shell synchrotron photons, so the 
IC emission contributes to 100 keV range as well as to the 10 GeV range. As a result, the 
photon indices for the Band function below a few MeV evolve to a ^ —1.25 and /3 ^ —2.35 
near the peak flux {a ^ —1.52 and f3 ~ —2.13 for the isotropic approximation) and return 
to the initial values later. Similarly, the spectrum in GRB 080916C becomes fiat in ef{e)- 



plot, when the GeV emission brightens (jAbdo et al. Il2009d ). Thus, an internal shock with 



an external Band/photospheric component appears generically able to explain the delayed 
onset of the GeV emission and the spectral evolution observed in Fermi LAT GRBs. 

While almost all the injected energy is released as photons in the former cases of Runl- 
Run4, the small e^j^j^ in Run5 leads to a slower energy release. Within the observer time 
< 6 s, 79% of the injected energy is released in the isotropic approximation. The integrated 
energy until the observer time 100(1 + z)Rq/T'^ /c ~ 100 s is 91% of the injected energy (we 
follow the photon emission from R = Rq to 100i?o)- The rest of the energy is dissipated 
mostly via adiabatic cooling. 



4. Summary of the numerical results 

The results of our time-dependent simulations for the internal dissipation model are 
summarized as follows. 

• The temporal evolution of the photon energy distribution inside the shell affects the 
resultant spectra seen by the observer. Especially for weak magnetic fields, the evo- 
lution of the IC component shows the limitation on the steady state approximation. 
Time-resolved spectra are necessary in order to verify expected features of the spectral 
evolution, such as the decrease of the cut-off energy due to 77-absorption, the growth 
of the IC component, and a spectral softening due to electron cooling. 

• Even in our one-zone approximation, the power-law spectrum is produced above the 77 
cut-off energy. The spectral evolution and photon escape during the electron injection 
are essential for this spectral shape. 



We confirm that the light curves follow the FRED shape (iFenimore et al.lll996l ). The 
pulse widths and peak times are affected by not only the curvature effect but also the 
intrinsic evolution of the energy distribution. The decrease of the cut-off energy leads 
to an earlier GeV peak time, and the gradual electron cooling causes the delayed onset 
and broad pulse profile of the optical emission. 

The possible deviation mechanisms from the simple Band function (other than IC 
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emission) are synchrotron emissions from secondary electron-positron pairs, electron 
heating via SSA, and late synchrotron emission after the injection ends. The latter 
two mechanisms yield a spectral bump in the low energy portion, while the secondary 
pairs enhance the synchrotron flux in the high energy range. 

The retarded growth of the IC component can cause the delayed onset of GeV emission. 
This may be a partial reason for the delayed onsets seen in the Fermi-LAT GRBs. 

The combination of the IC and late synchrotron components can produce the "extra 
spectral component" from keV to GeV as seen in the Fermi-LAT GRBs. This is a 
counter example to the claim that the IC emission cannot explain the extra components 
in the Fermi-LAT GRBs. 

When the IC and late synchrotron components are prominent, a spectral evolution 
from the Band function to a power-law dominated shape is seen. This is similar to the 
spectral evolution observed in GRB 090510. 

We confirm t he essentia l features of external photon models such as proposed by 
Toma et al. I (120091 120101 ) . The IC emission using external Band photons coming from 
inside regions can reproduce the delayed onset of GeV emission. The effect of the 
anisotropic emissivity adds an extra factor for the delay timescale. The evolution of 
the spectral indices is similar to those seen in GRB 080916C. The long-lasting tail of 
GeV emission in this model is also an interesting feature for the Fermi-LAT GRBs. 



Discussion 



Internal shocks, or more generally internal dissipation, have been consi dered as the 
main source to explain the GRB prompt emission f lPiran Il2005t iMeszaros Il2006[). but se veral 



problems, such as the narrow distribution of the peak energy gpeak (iPreece et al 



2000) etc., 



have led to consideration of photosphere models ( iPaczynski Ill986t iMeszaros fc Rees Il2000l ) 
in which most of the outflow energy is released as thermal photons from the photosphere 
More recently, dissipative photosphere m odels have been invoked to explain the spec t ra and 



high efficiency in the prompt emission ([Giannios fc Spruit 



20081 : iRvde et al. II2OIOI : lBeloborodovll2010l : IZhang et al. 



2011 



2007; 



loka at al. 



Vurm et al. 11201 ll ). However 



2007; Pe'er 



internal dissipation regions located outside the photosphere, including the classical internal 
shock model, are still attractive, as summarized in §H The emission from accelerated elec- 
trons in the dissipation region provides a natural mechanism for emitting Ge V photons as 
detected with Ferrni-LAT (as for the dissipative photosphere models, see, e.g.. iBeloborodov 
2OIOI : IVurm et al. II2OIII ). Also for the light curve shape, the photospheric radius may be 
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much smaller than ctvar^"^ (^var is the variability timescale), so that the photosphere model 
must attribute the broad pulse shape to the long-duration activity of the central engine. 
On the other hand, the internal dissipation mo dels can easily re produce the FRED shape 



(JKobayashi et al. 1119971 ) and the spectral lags fIShen et al.ll2005[ ). Although some fraction 



of GRBs such as G RB 090902B may be explained with a simple photospheric component 
(jZhang et al.ll201ll ). the internal dissipation in the outflows can play an important role not 
only inside the photosphere but also outside. 

In this paper we test two models for the extra-spectral components detected with Fermi- 
LAT; the first is the simple IC emission from accelerated electrons and the second is the 
up-scattering of external photons. Both models require an internal dissipation region in 
the outflow, where electron acceleration occurs. The possible delay time in the simple IC 
model is within the pulse width, so in this model one needs some additional reason to 
explain the ~ 4 s delay in GRB 080916C. However the external photon model has the 
additional advantage of being able to reproduce such delayed onset s of the GeV emission . 



We note that external shock mo dels with a high bulk-Lorentz factor (JGhisellini et al. 112010 



but as 



Maxham et al. 



Kumar fc Barniol Duran 1120101) can naturally explain the delayed onset of GeV emission. 



(120 111 ) have pointed out, the early-phase energy output estimated 
from the GBM data are not sufficient to explain the GeV fluxes in the prompt phase (see 
also, iHe et al. Il201ll ). The variability of the GeV emission needs to be tested with much 



high photon statistics in order to verify whether their origin is internal or external. 



The cascade emiss i ons initiated by hadronic interactions (e.g 



Bottcher fc Dermer 1998 



Gupta fc Zhang 1120071 : lAsano &: Inoue 1120071 : lAsano et al. Il2009bl ) are also one of the inter 
esting options for explaining the GeV photon emissio n, related to the potentia l role of GRB as 



sources of ultra high-energy cosmic rays (UHECRs) (IWaxman 



trino emission (jWaxman and Bahcall 1119971 : iGuetta et al. 



2001 



1995 



Vietri 



1995) and neu- 



Murase fc Nagataki II2006I ). 



In particular, the flat spectrum of the extra componen t in GRB 090902B is well explained 
by the hadronic model with a flducial parameter set fjAsano et al. Il2010l ). However, that 
GRB may be an exceptional case; the detection of GeV photons generally requires a large 
bulk Lorentz factor and large emission radius to avoid the internal 77-absorptio n . Such 



situations reduce the efficiency of the photomeson production. As lAsano et al. I (l2009al ) 



showed, a much larger proton luminosity than gamma-ray luminosity is required to produce 
the extra spectral component in GRB 090510. Although the high proton loading of the 
hadronic models is favorable for the GRB-UHECR scenario, the leptonic models considered 
in this paper are energetically less demanding and appear to provide a broader agreement 
with the observed spectra and time delays. The spectral evolution in the GeV-TeV range 
could provide essential clues for distinguishing the leptonic and hadronic models. This paper 
is thus both a test of purely leptonic GRB models based on SSC and EIC schemes, and also 
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a first step toward developing a time-dependent code to simulate emissions including the 
hadronic processes. We plan to compare the evolutions of the GeV-TeV spectrum for both 
the leptonic and hadronic models in near future. 

Concerning the classical internal shock model, in this paper we have not focused on the 
various problems mentioned with these, which have served as a motivation for photosphere 
models. We note that the minimum energy of electrons £^in' (conventionally expressed by 
the energy and number fractions of accelerated electrons) is treated as a free parameter in our 
simulations. Thus, our results do not address questions regarding the na rrow distribution 



of gp eak or the empirical ^peak-correlations such as the Yonetoku relation (JYonetoku et al 



20041 ). Similarly, the open problem of the low-energy spectral index a is not addressed here. 
In our results in §3.11 and 13.21 the shorter cooling time th an the dynamical ti mescale makes 



a 



-1.5, while the typical observed value is a = —1.0 fjPreece et al. Il2000l ). To solve the 



Qi-problem, several me c hanisms have been proposed, such as the Klein-Nishina effec t on SSC 
JPerishev et al. II2OO1I : ISosniak et al. II2OO9I: JNakar et al. I bood: IWang et al.lboOQh and the 



decay of magnetic fields (jPe'er fc Zhang II2006I ). as well as the effect of the su perposition of 



differe nt o ptical depth and te mperature regions in photospheric models, e.g., iBeloborodov 
(120 lOl ) and lVurm et al. I (120 111 ). Other issues related to the magnetic g eometry and expansion 
dynamics or to nuclear collisions effects (e.g.. lMeszaros &: Rees 11201 ll ) can affect the location 
of the photosphere or provide additional delay mechanisms. Also, in the dissip ation region, 
we c an expect magnetic turbulence due to t he Kelvin- Helmholtz instability ( Zhang et al. 



20091 ) or the Richtmyer-Meshkov instability (JMizuno et al. 1 1201 ll : llnoue et al 



201 ih . Such 



highly disturbed long-wave magnetic fields may interact with h igh energy electrons. Th e 
electron heating due to the turbulence may harden the spectrum (lAsano fc Terasawa II2009I ). 
Such effects expected in the dissipation region should be considered. 

The code discussed in this paper will be further developed to treat hadronic processes 
as well. In the near future, we plan to carry out simulations for various situations involving 
dissipative photospheres and internal or external dissipation or shock regions. In order to 
test with these a realistic picture of GRB, such calculations will need to adopt a wide range 
parameter sets for the mag netic field, bulk Lorentz factor and baryon loading (see e.g., 
loka Il2010l : lloka at al. II2OIII ). Moreover, the code will be useful for simulating emissions of 
other high-energy sources, such as active galactic nuclei, supernova remnant, and clusters of 
galaxies. 
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